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ABSTRACT 



Oh 



The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS (MUCHFUSS) aims at finding hot subdwarf 
stars with massive compact companions like massive white dwarfs (M > 1.0 M Q ), neutron stars or stellar mass black holes. The 
existence of such systems is predicted by binary evolution theory and recent discoveries indicate that they exist in our Galaxy. 
First results are presented for seven close binary sdBs with short orbital periods ranging from 0.21 d to 1.5 d. The atmospheric 
parameters of all objects are compatible with core helium-burning stars. The companions are most likely white dwarfs. In one case 
the companion could be shown to be a white dwarf by the absence of light-curve variations. However, in most cases late type main 
sequence stars cannot be firmly excluded. Comparing our small sample with the known population of close sdB binaries we show 
that our target selection method aiming at massive companions is efficient. The minimum companion masses of all binaries in our 
sample are high compared to the reference sample of known sdB binaries. 

£C) ' Key words, binaries: spectroscopic - stars: subdwarfs 



1. Introduction 

y—{ • Send offprint requests to: S. Geier, 

i-H ' e-mail: geier@sternwarte.uni-erlangen.de Subluminous B stars (sdBs) are core helium-burning stars with 

* Based on observations at the Paranal Observatory of the European very thin hydrogen envelopes and masses around O.5M (Heber 

t-H . Southern Observatory for programme number 081.D-0819. Based on 119861 see Heber 120091 for a review). A large fraction of the 

• • observations at the La Silla Observatory of the European Southern s dB stars (40 % to 8 %) are members of shor t period binaries 

; Observatory for programmes number 082.D-0649 and 084.D-0348. (Maxted et al. I200T1 Napiwotzki et al. I2004ab . Several studies 

. Based on observations collected at the Centra Astronomico Hispano were undertaken to determine the orbital parameters of subdwarf 

Aleman (CAHA) at Calar Alto, operated lointly by the Max-Planck , • j r j • j c rv nn » .u 1 rv j 

5— i T ... . _>. . ' ... ' \.. . , 3 . • J. , ... - binaries, and found periods ranging from 0.07 to more than 10 d 

Institut lur Astronomie and the Instituto de Astrohsica de Andalucia ., ' n fi > ^nm- « , 
■ ■ ' (CSIC). Based on observations with the William Herschel Telescope wlth a P eak atOlS to 1 .0 d (e.g. Edelmann et al. [20051 Morales- 
and the Isaac Newton Telescope operated both by the Isaac Newton Rueda et al - EQQ3). For close binary sdBs, common envelope 
Group at the Observatorio del Roque de los Muchachos of the (CE) ejection is the most probable formation channel. In this 
Instituto de Astrofisica de Canarias on the island of La Palma, Spain, scenario two main sequence stars of different masses evolve in 
Based on observations with the Southern Astrophysical Research a binary system. The heavier one will reach the red giant phase 
(SOAR) telescope operated by the U.S. National Optical Astronomy first and fill its Roche lobe. If the mass transfer to the companion 
Observatory (NOAO), the Ministrio da Ciencia e Tecnologia of the [ s dynamically unstable, a common envelope is formed. Due to 
Federal Republic of Brazil (MCT), the University of North Carolina friction the two stellar cores lose orbital energy, which is de- 
al Chapel Hill (UNC), and Michigan State University (MSU). Based ited wkhin ^ envel and leads t0 a short ening of the bi- 
on observations obtained at the Gemini Observatory, which is oper- 

ated by the Association of Universities for Research in Astronomy, 

Inc., under a cooperative agreement with the NSF on behalf of the Australian Research Council (Australia), Ministrio da Cilncia e 

Gemini partnership: the National Science Foundation (United States), Tecnologia (Brazil) and Ministerio de Ciencia, Tecnologa e Innovacin 

the Science and Technology Facilities Council (United Kingdom), Productiva (Argentina). This paper uses observations made at the South 

the National Research Council (Canada), CONICYT (Chile), the African Astronomical Observatory (SAAO). 
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Table 1. Solved binary systems. 



SDSS name short name 



other names 



SDSSJ002323. 
SDSS Jl 13840 
SDSSJ150513 
SDSSJ165404. 
SDSSJ172624. 
SDSSJ204613. 
SDSSJ225638. 



99-002953.2 
68-003531.7 
52+110836.6 
25+303701.7 
09+274419.3 
40-045418.7 
34+065651.0 



J0023-0029 
Jl 138-0035 
J1505+1108 
J1654+3037 
J 1726+2744 
J2046-0454 
J2256+0656 



PB5916 
PG 1136-003 
PG 1502+113 
PG 1652+307 
PG 1724+278 

PG 2254+067 



nary period. Eventually the common envelope is ejected and a 
close binary system is formed, which contains a core helium- 
burning sdB and a main sequence companion. If the companion 
has already evolved to a white dwarf (WD) when the red giant 
fills its Roche lobe, a close sdB+WD binary is formed (Han et 
al. 120021 [2003). Under certain conditions, two consecutive CE 
phases are possible as well. 

In general it is difficult to put constraints on the nature of the 
close companions to sdB stars. Since most of the binaries are 
single-lined, only lower limits have been derived from the bi- 
nary mass functions, which are in general compatible with main 
sequence stars of spectral type M or compact objects like white 
dwarfs. Only in special and hence rare cases can tighter con- 
straints be put on the nature of the companions. 

Subdwarf binaries with massive WD companions turned out 
to be candidates for supernova type la (SN la) progenitors be- 
cause these systems lose angular momentum due to the emission 
of gravitational waves and start mass transfer. This mass transfer, 
either from accretion of He onto the WD during the sdB phase 
(e.g. Yoon & Langer 120041 and references therein), or the sub- 
sequent merger of the system after the sdB star itself has turned 
into a WD (Tutukov & Yungelson [T98TJ Webbink [19841 may 
cause the companion to approach the Chandrasekhar limit and 
explode as SN la. 

SN la play a key role in the study of cosmic evolution (e.g. 
Riess et al. [19981 Leibundgut I200T] Perlmutter et al. \T§9% . 
One of the best known candidate systems for the double degen- 
erate merger scenario is the sdB + WD binary KPD 1930+2752 
(Maxted et al. 2000a] Geier et al. l2007b . Mereghetti et al. ( [20091 
showed that in the X-ray binary HD 49798 a massive (> 1 .2 M ) 
white dwarf accretes matter from a closely orbiting subdwarf O 
companion. The predicted amount of accreted material is suffi- 
cient for the WD to reach the Chandrasekhar limit. This makes 
HD 49798 another candidate for SNIa progenitor. Furthermore, 
Perets et al. ( 120 101 ) showed that helium accretion onto a white 
dwarf may be responsible for a subclass of faint and calcium- 
rich SN lb events. 

Geier et al. d2008l I2010al 1201 Obi ) analysed high resolution 
spectra of sdB stars in close binaries. Assuming synchronised 
rotation they constrained the masses and the nature of the unseen 
companions in 31 cases. While most of the derived companion 
masses were consistent with either late type main sequence stars 
or white dwarfs, the compact companions of some sdBs may be 
either massive white dwarfs, neutron stars (NS) or stellar mass 
black holes (BH). However, Geier et al. (2010b) also showed that 
the assumption of orbital synchronisation in close sdB binaries 
is not always justified and that their analysis suffers from huge 
selection effects. 

The existence of sdB+NS/BH systems is predicted by bi- 
nary evolution theory (Podsiadlowski et al. 120021 Pfahl et al. 



2003). The formation channel includes two phases of unsta- 
ble mass transfer and one supernova explosion. The fraction of 
sdB+NS/BH systems is predicted to be about 2% of the close 
sdB binaries (G eier et al. l2010bl >. Yungelson & Tutukov ([20531 ) 
and Nelemans (2010) performed independent binary evolution 
calculations and confirm that sdB+NS/BH systems should ex- 
ist. According to the results of Nelemans (2010) about 1% of 
the subdwarfs in close binaries should have a neutron star com- 
panion, whereas only 0.01% should be orbited by a black hole. 
Yungelson & Tutukov (2005) predict the sdB+NS fraction to be 
of the order of 0.8%. 

Since sdB stars eventually evolve to WDs there should also 
exist a population of white dwarfs with massive compact com- 
panions. Badenes et al. (2009) reported the discovery of a close 
binary consisting of a massive white dwarf and an unseen neu- 
tron star or black hole companion, but Marsh et al. (2010) most 
recently showed that the system is double-lined and consists of a 
massive white dwarf orbited by a low mass white dwarf. The sys- 
tem mass is below the Chandrasekhar limit. Their results were 
confirmed by Kulkarni & van Kerkwijk (120101 ). Common enve- 
lope ejection was proposed as the most likely formation channel 
for the binary PSR J1802-2124, which consists of a millisecond 
pulsar and a CO white dwarf in close orbit (P = 0.7 d, Ferdman 
et al. 120101 1. This peculiar system may have evolved through an 
earlier sdB+NS phase. 



2. The MUCHFUSS project 

The discovery of sdB binary candidates with massive compact 
companions provides a first hint that a whole population of non- 
interacting binaries with such companions may be present in our 
Galaxy. The known candidate sdB+NS/BH binaries have low or- 
bital inclinations (15 - 30°, Geier et al. 2010b). High inclination 
systems must exist as well and should be more numerous. In this 
case a determination of the orbital parameters is sufficient to put 
a lower limit to the companion mass by calculating the binary 
mass function. If this lower limit exceeds the Chandrasekhar 
mass and no sign of a companion is visible in the spectra, the 
existence of a massive compact companion is proven without 
the need for any additional assumptions. 

The project Massive Unseen Companions to Hot Faint 
Underluminous Stars from SDSSQ (MUCHFUSS) aims at find- 
ing sdBs with compact companions like massive white dwarfs 
(M > 1.0 M Q ), neutron stars or black holes. About 80 binaries 
have been selected for follow-up. Survey and target selection are 
described in detail in Geier et al. (2010c). The same selection cri- 
teria that we applied to find such binaries are also well suited to 
single out hot subdwarf stars with constant high radial velocities 
(RV) in the Galactic halo and search for hypervelocity stars. First 
results of this second part of the project (Hyper-MUCHFUSS) 
are presented in Tillich et al. (1201 Oi l. 

Here we present the spectroscopic analysis of the first sdB 
binaries discovered in the course of the MUCHFUSS project 
(see Table Q]). In Sect. [3] the observations and the data reduc- 
tion are described. Sects. [4]and[5]deal with the determination of 
the orbital and atmospheric parameters of the sdB stars. Sect. [6] 
explains the way the minimum masses of the unseen compan- 
ions are constrained, while results are presented in Sect. [7] The 
efficiency of our target selection is discussed in Sect. [8] a short 
summary and an outlook are eventually given in Sect. [9] 



Sloan Digital Sky Survey 



2 



S. Geier et al.: Massive unseen companions to hot faint underluminous stars from SDSS (MUCHFUSS) 



(J 
+ 




J0023-0029 (EF0SC2) 




J0023-0029 (SDSS) 






J1 138-0035 (IDS) 



J1 505+1 108 (TWIN) 




J 1654 + 3037 (SDSS) 



J1 726 + 2744 (ISIS) 
J2046-0454 (F0RS1) 



J2256+0656 (ISIS) 




_L 



_i i i_ 



_i i i_ 



4000 



4500 

x[A] 



5000 



Fig. 1. Medium resolution spectra of the programme stars taken with different instruments. Multiple observations of the same target 
have been shifted to rest wavelength and coadded. 



3. Multi-site observations and data reduction 

Follow-up medium resolution spectra were taken during de- 
dicated follow-up runs (see Table [2} with the EFOSC2 spec- 
trograph (R =i 2200, A = 4450 - 5 110 A) mounted at the 
ESO NTT, the ISIS spectrograph (R 4000, A = 3440-5270 A) 
mounted at the WHT, the TWIN spectrograph mounted at the 
CAHA-3.5m telescope (R » 4000,^ = 3460 - 5630 A), the 
Goodman spectrograph mounted at the SOAR telescope (R 
2500,/! = 3500 - 6160 A), the GMOS spectrograph (R 
1200,/! = 3770 - 4240 A) mounted at the Gemini North tele- 
scope and the IDS spectrograph mounted at the Isaac Newton 
Telescope (R a 1400,^ = 3000 - 6800 A). Informations about 
data taken in the course of our survey are provided in Geier et 



al. (1201 Pel ). Additional data could be gathered, when our tar- 
gets were observed with the IDS spectrograph (March 2007, ob- 
server: T. M., C. C; R 4000, /I = 3930 - 5100 A) and the 
grating spectrograph (March 2003, April 2004, observer: T. M.; 
R - 4600, A = 4170 - 5030 A) mounted at the 1.9m Radcliffe 
Telescope. Example spectra are shown in Fig.[TJ 

In order to obtain a good wavelength calibration, arc lamp 
exposures have been taken before or after the single exposures. 
In addition to that bright single sdBs have been taken as RV stan- 
dards in most of the runs. In some cases the RVs of certain instru- 
ments (TWIN, GMOS) had to be corrected by a constant offset 
of up to =i 50kms _1 , which was derived from the RV measure- 
ments of the standard stars. The slit width was always chosen 
to be smaller than the size of the seeing discs to minimize sys- 
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Fig. 2. Radial velocity plotted against orbital phase. The RV data were phase folded with the most likely orbital periods. The 
residuals are plotted below. The RVs were measured from spectra obtained with SDSS (rectangles), CAHA3.5m/TWIN (upward 
triangles), WHT/ISIS (diamonds), INT/IDS (downward triangles), ESO-VLT/FORS 1 (triangles turned to the left), Gemini/GMOS 
(triangles turned to the right), ESO-NTT/EFOSC2 (circles), SOAR/Goodman (hexagons) and SAAO-1.9m/Grating (stars). 

Table 3. Derived orbital parameters. 



Object 



To 

[-2450000] 



P 

[d] 



[kms- 1 ] 



K 

[kms" 1 ] 



J0023-0029 
Jl 138-0035 
J1505+1108 
J1654+3037 
J1726+2744 
J2046-0454 
J2256+0656 



5069.850 ± 0.008 
4991.388 +0.001 
4938.867 ± 0.002 
4991.5322 ±0.0008 
4981.667 ±0.005 
4693.352 ±0.002 
5070.662 ± 0.002 



1.4876 ±0.0001 
0.207536 ± 0.000002 
0.74773 ± 0.00005 
0.25357 ± 0.00001 
0.50198 ± 0.00005 
0.24311 ±0.00001 
0.7004 ± 0.0001 



16.4 ±2.1 
23.3 ± 3.7 

-77.1 ± 1.2 

40.5 ± 2.2 
-36.7 ±4.8 

87.6 ±5.7 
-7.3 ±2.1 



81.8 ±2.9 

162.0 ± 3.8 
97.2 ± 1.8 

126.1 ±2.6 
118.9 ± 3.7 
134.3 ±7.8 
105.3 ±3.4 
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Fig. 3. Radial velocity curves (see Fig[2j. 



tematic errors due to movement of the objects within the slit. 
Reduction was done either with the MIDAS, IRAF or PAMELifl 
and MOLLY 2 packages. 



4. Orbital parameters 

The radial velocities were measured by fitting a set of mathemat- 
ical functions (Gaussians, Lorentzians and polynomials) to the 
hydrogen B aimer lines as well as helium lines if present using 
the FITSB2 routine (Napiwotzki et al. I2004bb . The RVs of the 
GMOS spectra have been measured by fitting three Gaussians 
to the H y line. Three functions are used to match the contin- 
uum, the line and the line core, respectively and mimic the typ- 



2 http://www2.warwick.ac.uk/fac/sci/physics/research/astro/people 
/marsh/software 



ical Voigt profile of spectral lines. The profiles are fitted to all 
suitable lines simultaneously using ^-minimization and the RV 
shift with respect to the rest wavelengths is measured. The RVs 
and formal lcr-errors are given in Appendix IB! Assuming cir- 
cular orbits sine curves were fitted to the RV data points in fine 
steps over a range of test periods. For each period the x 2 of the 
best fitting sine curve was determined. The result is similar to 
a power spectrum with the lowest x 1 indicating the most likely 
period (see Fig.|4j». In order to estimate the significance of the or- 
bital solutions and the contributions of systematic effects to the 
error budget, we normalised the x 1 of the most probable solution 
by adding systematic errors in quadrature until the reduced x 2 
reached ^ 1 .0. Using these modified uncertainties we performed 
Monte Carlo simulations for the most likely periods. For each 
simulation a randomised set of RVs was drawn from Gaussian 
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Table 4. Significance of the circular orbital solutions. The best solutions for the orbital periods are given together with their minimum 
X 1 and reduced^ 2 values as well as the number n of RVs. The second best aliases (further than 1 % away from the best solution) and 
the A^ 2 -values with respect to the best solutions are given as well. The systematic error adopted to normalise the reduced^ 2 (e nor m) 
is given for each case. The probabilities for the orbital period to deviate from our best solution by more than 1% (f>f a ise[l%]) or 10% 
(ffaise[10%]) are given in the last columns. 



Object 


best solution 


x 2 


reduced 


2nd best alias 


A* 2 


n 


^-norm 


log Pfalse[l%] 


logPfalse[10%] 




[d] 






[d] 






[km s- 1 ] 






J0023-0029 


1.4876 


157 


3.74 


0.5976 


130 


47 


8.0 


-3.0 


-3.4 


Jl 138-0035 


0.207536 


213 


5.33 


0.260192 


426 


45 


16.0 


-3.5 


-3.5 


J1505+1108 


0.74773 


155 


4.30 


0.75709 


679 


41 


7.0 


< -4.0 


< -4.0 


J1654+3037 


0.25357 


18 


0.54 


0.20397 


64 


38 




< -4.0 




J 1726+2744 


0.50198 


82 


2.48 


1.00998 


77 


38 


12.0 


-1.2 


-1.9 


J2046-0454 


0.24311 


52 


3.05 


0.31971 


39 


22 


17.0 


-1.1 


-1.1 


J2256+0656 


0.7004 


276 


6.13 


2.1903 


976 


50 


13.0 


< -4.0 


< -4.0 



Table 2. Follow-up observations 2009/2010. The first column 
lists the date of observation, while in the second the used tele- 
scope and instrumentation is shown. In the third column the ob- 
servers are listed. 



Date 



Telescope & Instrument Observer 



2009/06/05-2009/06/09 ING-INT/IDS 



2009/07/22-2009/07/26 
2009/08/24-2009/08/27 
2009/11/08-2009/11/12 
April/August 2009 
2010/02/12-2010/02/15 



CAHA-3.5m/TWIN 

ING-WHT/ISIS 

ESO-NTT/EFOSC2 

Gemini-North/GMOS 

SOAR/Goodman 



R. 
T. O. 
T. K. 
S. G. 
T. K. 
Service 
B. B. 



R. O.. 



distributions with central value and width corresponding to the 
RV measurements and the analysis repeated. From these simula- 
tions the probabilities for the orbital periods to deviate from our 
best solution by more than 1% or 10% were calculated. 

In order to derive most conservative errors for the RV semi- 
amplitude K and the system velocity y we fixed the most likely 
period and created new RV datasets with a bootstrapping algo- 
rithm. Ten thousand RV datasets were obtained by random sam- 
pling with replacement from the original dataset. In each case an 
orbital solution was calculated in the way described above. The 
standard deviation of these results was adopted as error estimate. 
The RV curves are given in Figs.[2]and[3] The residuals of the RV 
curves after subtracting the best orbital solution are of the same 
order in all cases (see Figs. 12 El). The accuracy is limited by the 
resolution of the spectra and their signal-to-noise. Combining 
data obtained with different instruments is also expected to con- 
tribute to the systematic error. Nevertheless, we found that all 
orbital solutions given here are significant (see Table|4]l. 

Edelmann et al. (120051) reported the discovery of small ec- 
centricities (<? < 0.06) in the orbital solutions of five close hot 
subdwarf binaries. All of these binaries are expected to have 
formed via common envelope ejection. Although the CE phase 
is very short, it should nevertheless be very efficient in circular- 
ising the binary orbits. That is why the discovery of Edelmann et 
al. ((2005) came as a surprise. Napiwotzki et al. (in prep.) found 
more such systems with even shorter periods. 

In order to investigate whether the orbital solutions of our 
programme binaries can be improved by allowing for eccen- 
tricity, we fitted eccentric orbits to our radial velocity data and 



performed statistical tests (F-test, see Pringle 119751 and the 
Bayesian information criterion BIC) to check whether eccentric 
solutions are significant or not. In all cases the circular solutions 
were preferred. However, the derived upper limits for the orbital 
eccentricities range from 0.15 to 0.3, which means that low ec- 
centricities as the ones reported by Edelmann et al. (2005) can- 
not be firmly excluded. 

5. Atmospheric parameters 

Atmospheric parameters have been determined by fitting model 
spectra to the hydrogen Balmer and helium lines in the way 
described in Geier et al. (120071) . The single spectra have been 
corrected for their orbital motion and coadded. Depending on 
the effective temperature of the stars, LTE models with so- 
lar metallicity (r e ff < 30 000 K) or ten times solar metallicity 
(T e g > 30 000 K) have been used. The enhanced metallicity 
models account for the radiative levitation of heavy elements in 
the diffusion dominated atmospheres (for a detailed discussion 
see O'Toole & Heber|2006]). 

In order to investigate systematic effects introduced by the 
individual instruments, especially the different resolutions and 
wavelength coverages, the parameters have been derived sepa- 
rately from spectra taken with different instruments. As can be 
seen in Table I A. 11 no constant systematic shifts are present. The 
weighted means have been calculated and adopted as final so- 
lutions. Typical systematic errors introduced by different model 
grids are of the order of +0.05 in logg and 500 K in T e g (e.g. 
Lisker et al. 120051 Geier et al. 120071 ) . These uncertainties were 
added in quadrature to the statistical errors. 

Three of our programme stars have been classified as hot 
subdwarfs by Eisenstein et al. (2006), but the authors pointed 
out that the atmospheric parameters of the sdO/Bs given in their 
catalogue are not accurate. 

All stars of our sample are situated on or near the Extreme 
Horizontal Branch (EHB) and are most likely core-helium burn- 
ing stars (see Fig. [5]). Since the orbital periods of these binaries 
are short, they can only have formed via common envelope ejec- 
tion. Population synthesis models (Han et al. 2002, 2003) pre- 
dict a mass range of M^b = 0.37 - 0.48 M for sdBs in binaries 
formed in this way. The mass distribution shows a sharp peak 
at a mass of about 0.47 M . This theoretical mass distribution is 
consistent with analyses of close binary systems (e.g. Geier et al. 
120071 For et al. 120101 ) as well as asteroseismic analyses of pul- 
sating sdBs (see Charpinet et al. 2008 and references therein). If 
the progenitor star was massive enough on the main sequence to 
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Fig. 4. x 2 plotted against orbital period. The lowest peak corresponds to the most likely solution. 



ignite core helium-burning under non-degenerate conditions, the 6. Constraining the nature Of the unseen 

sdB mass may be as low as 0.3 M Q . A small fraction of the sdB companions 
population is predicted to be formed in that way (Han et al. 120021 

2003 ). Especially for sdB binaries with massive companions this since me programme stars are single-lined spectroscopic bina- 



formation scenario may become important. 



ries, only their mass functions can be calculated. 



fm 



(Mcomp + M sdB ) 2 2ttG 



(1) 



Although the RV semi-amplitude K and the period P can be 
derived from the RV curve, the sdB mass M s( jb, the companion 
mass M comp and the inclination angle i remain free parameters. 
Adopting M s dB = 0.47 M and i < 90° we derive a lower limit 
for the companion mass (see Table|5]l. 
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Fig. 5. T e ff - log ^-diagram. The helium main sequence (HeMS) 
and the EHB band (limited by the zero-age EHB, ZAEHB, and 
the terminal-age EHB, TAEHB) are superimposed with EHB 
evolutionary tracks from Dorman et al. (1993 ). 

For minimum companion masses lower than 0.45 M the 
companion may be a late type main sequence star or a compact 
object like a WD. Main sequence stars in this mass range are 
outshined by the sdBs and not visible in optical spectra (Lisker 
et al. I20051 >. That is the reason why the companions' nature still 
remains unknown for most of the ^80 known sdB systems with 
low minimum companion masses (see Fig. |7|. If on the other 
hand the minimum companion mass exceeds 0.45 M , spectral 
features of a main sequence companion become visible in the 
optical. The non-detection of such features therefore allows us 
to exclude a main sequence star. The companion must then be 
a compact object. More massive compact companions like mas- 
sive WDs, neutron stars or black holes are more likely as soon as 
the minimum mass exceeds 1 .00 M or even the Chandrasekhar 
limit 1. 40 M . 

Due to the fact that we selected targets with high RV shifts, 
the distribution of orbital inclinations in our target sample is 
not random any more. Our selection strategy strongly favours 
high inclination angles, and therefore the companion masses are 
likely to be close to their minimum values. The probability of 
detecting eclipses, reflection effects or variations caused by el- 
lipsoidal deformation in the light curves of systems with short 
orbital periods should therefore be significantly higher than in 
an unbiased sample. 

7. Results 

The spectra of all stars in our sample have been checked for 
spectral features of their companions. Hot subdwarfs with faint 
main sequence companions usually show spectral lines of the 
Mg i triplet at a 5170 A (Lisker et al. 120051 1 and the Can triplet 
at a 8650 A. No such features are visible in the spectra of our 
programme stars (see e.g. Fig. [TJ. Stark & Wade (2003) ana- 
lysed optical and IR photometry (2MASS) and found no indi- 



Table 5. Derived minimum masses and most probable nature of 
the companions. 



Object 


f(M) 


Af 2rai „ 


Companion 




[M ] 


[M Q ] 




J0023-0029 


0.084 


0.40 


MS/WD 


Jl 138-0035 


0.091 


0.42 


WD 


J1505+1108 


0.071 


0.37 


MS/WD 


J1654+3037 


0.053 


0.32 


MS/WD 


J 1726+2744 


0.087 


0.41 


MS/WD 


J2046-0454 


0.061 


0.34 


MS/WD 


J2256+0656 


0.085 


0.40 


MS/WD 



cation of an IR-excess caused by a cool companion in the case 
of J1654+3037. According to the catalogue of Reed & Stiening 
(2004), who performed a similar analysis, J1505+1108 shows 
signs of an IR-excess in the H and K-bands, but the large errors 
of these measurements and the missing spectral signatures of a 
cool companion in the SDSS spectra are strong indications, that 
no visible companion is present. 

J1654+3037 and J2046-0454 have very similar orbital pa- 
rameters. The periods are short (0.25 d) and the minimum com- 
panion masses are constrained to 0.32 M and 0.34 M . Whether 
the companions are M dwarfs or WDs is therefore not yet clear. 
In the former case a reflection effect should be easily detectable 
in the light curves. Photometric follow-up will allow us to clarify 
the nature of the companions. 

The companion of the short period (0.2 d) system 
Jl 138-0035 is most likely a white dwarf. The minimum com- 
panion mass is constrained to 0.42 M and no sign of a compan- 
ion is seen in the spectra. A light curve taken by the SuperWASP 
project (Pollacco et al. 2006) shows no variation exceeding 
a 1% (see Fig. [6]). Due to the short period of this system a reflec- 
tion effect should be visible, if the companion should be a cool 
main sequence star. The absence of such a variation leads to the 
conclusion that the companion is most likely a white dwarf. 

The orbital periods of J 1726+2744 (0.5 d), J2256+0656 
(0.7 d) and J1505+1108 (0.74d) are longer. Their minimum 
companion masses are similar (0.37 - 0.41 M ) and close to 
the border between main sequence stars and white dwarfs. The 
companions of J1726+2744 and J2256+0656 are most likely 
WDs. Koen ( 120091 ) and Shimanskii et al. ( 120081 ) recently showed 
that reflection effects can still be detected in the light curves of 
sdB binaries with similar orbital periods. A reflection effect in 
J0023-0029 on the other hand is most likely not detectable, be- 
cause the orbital period is too long (1 .5 d). 

8. Efficiency of target selection 

The goal of the MUCHFUSS project is to find sdB binaries with 
massive compact companions and study this population of close 
binaries. We tried to optimise our target selection to achieve 
this goal. Fig. [7] illustrates the efficiency of our target selection. 
The RV semiamplitudes of all known sdB binaries with spec- 
troscopic solutions (open symbols) are plotted against their or- 
bital periods (Geier et al. 1201 Pel) . Binaries which have initially 
been discovered in photometric surveys due to indicative fea- 
tures in their light curves (eclipses, reflection effects, ellipsoidal 
variations) are marked with open circles. Binaries discovered by 
RV variations from time resolved spectroscopy are marked with 
open diamonds. The dashed, dotted and solid lines mark the re- 
gions to the right where the minimum companion masses derived 
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0.1 - 




-0.05 - 



-0.1 - 



-0.15 I ' 1 ' 1 ' ' 

-10 12 
phose 

Fig. 6. SuperWASP light curve of Jl 138-0035 folded to the or- 
bital phase. The 11213 data points taken between 2006/07/05 
and 2009/07/02 are binned to 100 phase bins. Relative flux is 
plotted against the orbital phase. 



from the binary mass function (assuming 0.47 M Q for the sdBs) 
exceed 0.45 M , 1 .00 M and 1 .40 M . 

Most of the known sdB binaries are situated beneath the 
0.45 M line, which means that the companion type cannot be 
constrained from the mass function alone. Photometry is neces- 
sary to clarify the companions' nature in these cases. The most 
massive sdB binary known to date is KPD 1930+2752 with a 
WD companion of 0.9 M . This short period system has been 
discovered based on indicative features in its light curve (upper 
left corner in Fig. [7j Billeres et al. 2000). 

The seven binaries from the MUCHFUSS project are marked 
with filled diamonds. It can be clearly seen that they belong 
to the sdB binary population with the largest minimum masses 
close to 0.45 M . We therefore conclude that our target selection 
is efficient and singles out sdB binaries with massive compan- 
ions. 



9. Summary and Outlook 

A multi-site follow-up campaign is being conducted with 
medium resolution spectrographs mounted at several different 
telescopes of mostly 2 m to 4m-class. First results were pre- 
sented for seven close binary sdBs with short orbital periods 
ranging from ^ 0.21 d to 1.5 d and most likely compact com- 
panions. The atmospheric parameters of all objects are compati- 
ble with core helium-burning stars on the EHB. Comparing our 
small sample with the known population of close sdB binaries 
we are able to show that our target selection method is efficient. 
All binaries solved up to now have high minimum companion 
masses compared to the rest of the sdB binary population. 

Up to now we have found significant orbital solutions for 
about 10% of our target sample. Photometric follow-up obser- 
vations will allow us to clarify the nature of the companions in 
most cases. A database of more than 700 spectra has been built 




o I i i 1 1 1 i i i I i i i I i i i 

0.1 1 10 100 

p [d] 



Fig. 7. The RV semiamplitudes of all known sdB binaries 
with spectroscopic solutions plotted against their orbital periods 
(Geier et al. l2010cl >. Binaries which have initially been discov- 
ered in photometric surveys due to indicative features in their 
light curves (eclipses, reflection effects, ellipsoidal variations) 
are marked with open circles. Binaries discovered by detection 
of RV variations from time resolved spectroscopy are marked 
with open diamonds. The dashed, dotted and solid lines mark 
the regions to the right where the minimum companion masses 
derived from the binary mass function (assuming 0.47 M for the 
sdBs) exceed 0.45 M , 1 .00 M and 1 .40 M . The seven binaries 
from the MUCHFUSS project are marked with filled diamonds. 



up and some binaries will be solvable with only a few additional 
RV points. 
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Appendix B: Radial velocities 

Table B.l. J0023-0029 



mid-HJD 


RV [kms- 1 ] 


Instrument 


-2450000 






1816.782390 


92.0 ± 10.0 


SDSS 




94.1 ± 6.0 




1816.806476 


92.7 ± 7.0 




1885.575567 


5.6 ± 9.0 




1885.587847 


2.2 ± 11.0 




1885.600718 


2.4 ± 15.0 




1885.614294 


-3.1 +9.0 




1885.627153 


7.7 ± 9.0 




1885.639444 


-9.0 ± 8.0 




1899.578750 


-36.7 ± 14.0 




1899.590972 


-17.2 ±9.0 




1899.603316 


-27.7 ± 10.0 




1899.615509 


-12.1 ± 11.0 




1900.573652 


-13.8 + 7.0 




1900.585712 


-9.6 ± 22.0 




1900.599109 


-25.3 + 16.0 






90.6 + 3.0 


ISIS 


5068.62614 


94.8 + 3.7 




5068.63335 


O 1 A 1 A f\ 

ol.4 ± 4.1) 




5068.65013 


1 eft r\ 1 a r\ 

102.0 + 4.0 




5068.65736 


1 m C 1 A A 

103.5 + 4.4 




5068.66457 


94.9 ± 3.7 




5069.69141 


-25.4 + 4.6 




5069.69862 


-65.8 ±4.1 




5069.70585 


-32.1 ± 2.9 




5069.71578 


-28.2 ± 2.6 




5069.72300 


-17.9 ±4.3 




5069.73023 


-6.2 ±4.1 




5071.59347 


71.9 ±5.9 




5071.60068 


79.3 ± 6.4 




5071.60442 


84.1 ±4.2 




5071.61757 


91.7 ±5.8 




5071.62478 


100.7 ± 2.6 




5071.63199 


93.9 ± 3.2 




5144.55811 


101.0 ± 10.1 


EFOSC2 


5144.59234 


97.1 ± 10.1 




5144.63989 


76.5 ± 12.6 




5145.52594 


-28.5 ± 9.9 




5145.61524 


-16.9 ±7.8 




5145.65247 


-7.1 ± 10.2 




5146.54220 


2.5 ± 8.0 




5146.64161 


-27.0 ±8.1 




5146.69219 


-57.2 ± 8.9 




5147.59722 


82.4 ± 10.4 




5147.64906 


81.5 ±9.5 




5147.66103 


85.2 ±7.2 




5147.67897 


80.0 ± 9.8 





Table B.2.J1 138-0035 



mirl-HIT") 

1111 Q LIJLJ 


RV [kms"'] 




-2450 000 








1629.831447 


-16.0 ± 


33.7 


SDSS 


1629.861748 


-131.2 ±5.8 




1 679 875990 


-139.6 ± 


10.5 




1 679 8901 97 


-124.6 ± 


11.9 




1 fi^O 84QS4Q 

1 \J J\J. Ot^7J^7 


42.3 ± 


12.3 
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-8.6 ± 5.8 




1658.666991 


16.0 ±6.1 




1658.723102 


-127.9 ± 2.4 
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-120.4 : 


t7.4 
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-73.0 ± 


32.8 


SAAO 
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IDS 
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-76.5 : 
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^77 1 .T'JJZ.J 
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^77 1 .TJUO / 
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477 1 .^f J 1 jl 
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Goodman 


5240.64678 
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Table B.3. J 1505+ 1108 



mid-HJD 


RV [kms '] 


Instrument 


1 A CA f\f\f\ 

-2450000 






3848.858414 


5.3 ± 12.0 


SDSS 


3848.906794 


27.7 ± 12.2 




3849.863669 


-126.6 ± 10.0 




3850.893113 


-96.1 ± 8.1 






-4o.U ± o.U 


TWTM 


4600.487332 


-65.4 ± 10.0 




4692.348022 


22.2 ± 13.0 




4694.404665 


-44.7 ± 9.7 




4696.399472 


-148.2 ± 10.8 




4980.44927 


-130.0 ± 9.2 




4980.53561 


-175.5 ± 9.9 




4981.57813 


-34.5 ± 13.8 




4982.55136 


5.4 ± 12.1 




4936.672663 


-37.3 ± 2.8 


GMOS 


4936.676403 


-31.0 ±3.0 




4936.680142 


-34.7 ± 3.0 




4936.683881 


-32.5 ±3.1 




4936.871029 


-3.5 ±2.8 




4936.874768 


-2.1 ± 2.8 




4936.878512 


0.9 ± 2.8 




4936.882253 


-5.5 ± 2.9 




4937.64103 


0.8 ±3.1 




4937.64477 


7.0 ± 2.9 




4937.64852 


6.1 ±3.0 




4937.65226 


4.4 ± 2.9 




4937.85083 


-158.7 ±2.9 




4937.85457 


-161.9 ±2.9 




4937.85831 


-158.5 ±2.9 




4937.86206 


-158.7 ± 3.0 




4938.75627 


-150.6 ±2.5 




4938.76001 


-149.4 ± 2.5 




4938.76375 


-156.3 ±2.5 




4938.76749 


-146.0 ± 2.5 




4939.67990 


-31.1 ±2.8 




4939.68364 


-24.8 ± 2.9 




4939.68738 


-26.7 ± 2.8 




4939.69112 


-23.6 ± 2.8 




4943.63233 


-2.9 ± 3.0 




4943.63607 


-9.2 ± 2.9 




4943.63981 


-7.8 ± 2.9 




4943.64355 


-10.3 ± 2.7 





Table B.4. J1654+3037 



mid-HJD 


RV [km s 1 ] 


Instrument 


-2450 000 






2789.917095 


119.4 ± 9.7 


SDSS 


2789.933032 


84.5 ± 7.7 




2790.913235 


146.7 ± 7.2 




2790.929502 


130.3 ± 8.5 




4586.567656 


-6.0 ±7.8 


TWIN 


40VZ.;5u / J ly 


IJJ.Z ± o.U 




4693.380826 


148.1 ± 8.1 




4694.433521 


155.9 ± 8.0 




5037.44285 


70.1 ± 10.6 




5037.47019 


150.9 ± 8.4 




5037.50213 


156.2 ± 11.1 




5038.42271 


-21.8 ± 11.0 




5038.48616 


155.0 ± 10.8 




5038.49857 


161.0 ± 7.9 




4988.47623 


-4.3 ± 25.5 


IDS 


4988.49036 


56.8 ± 16.6 




4988.50437 


97.9 ± 17.5 




4988.52125 


129.9 ± 15.6 




4988.53530 


172.6 ± 17.3 




4988.54942 


179.7 ± 18.0 




4988.56430 


150.0 ± 16.0 




4991.47179 


-84.8 ± 17.7 




4991.47888 


-82.0 ± 17.1 




4991.48596 


-92.2 ± 16.6 




4991.49305 


-58.5 ±7.1 




4991.50013 


-62.6 ± 17.5 




4991.50927 


-44.1 ± 16.5 




t?7 1 .J 1 / JO 


1 1 9 + 1 S 1 

-11. Z. IE U.l 




4991.52577 


23.0 ± 18.4 




4991.53402 


33.3 ± 16.2 




4991.54224 


66.3 ± 16.0 




4991.55155 


103.8 ± 30.1 




4991.55981 


128.0 ± 17.6 




4991.56805 


142.9 ± 16.5 




4991.57629 


161.2 ± 17.1 




4991.58453 


178.5 ± 16.4 
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Table B.5. J1726+2744 



mid-HJD 


RV [kms '] 


Instrument 


r\ a c f\ r\r\r\ 

-2450000 






3905.819525 


-168.9 ± 10.9 


SDSS 


3905.833513 


-157.3 ± 10.7 




3905.853553 


-133.2 ± 10.5 




3905.866007 


-123.5 ± 15.8 




4979.59252 


-125.2 ± 14.4 


TWIN 


4979.63382 


-98.5 ±11.6 




4980.47951 


-155.4 ± 12.9 




4980.58604 


-141.2 ± 11.4 




4981.53132 


-120.0 ± 8.0 




4981.58892 


-112.5 ± 12.9 




4981.62939 


-46.7 ± 12.0 




4981.64123 


-53.1 ±9.7 




4981.65164 


-35.2 ± 16.9 




4982.56453 


-154.2 ± 11.2 




4982.65151 


-65.0 ±21.6 




4983.53806 


-152.7 ± 11.5 




4983.54720 


-176.7 ± 12.0 




4983.55638 


-150.5 ± 15.2 




4983.56557 


-146.3 ± 24.3 




/1QQQ ^a/iAQ 

4Voj.j /4oo 






5037.45679 


39.0 ± 12.0 




5037.48545 


62.0 ± 12.7 




5038.43771 


38.0 ± 16.1 




5038.45381 


58.0 ± 12.9 




5038.47045 


80.0 ± 13.0 




5039.49762 


69.0 ± 12.7 




5039.52128 


87.0 ± 12.0 




5039.53838 


85.0 ± 12.9 




4992.48445 


-99.5 ± 23.1 


IDS 


4992.49617 


-114.1 ± 18.5 




4992.50788 


-81.2 ± 13.7 




5068.47186 


-92.7 ± 5.2 


ISIS 


5068.47906 


-79.4 ± 4.4 




5068.48628 


-73.9 ± 5.3 




5069.45355 


-120.3 ± 3.2 




5069.46426 


-109.9 ± 2.7 




5069.46786 


-98.2 ± 4.5 




5069.47855 


-79.9 ± 3.8 





Table B.6. J2046-0454 



mid-HJD 


RV[kms-'] 


Instrument 


-2450 000 






3269.661429 


109.6 ± 13.3 


SDSS 


3269.675556 


128.1 ± 9.6 




3269.691435 


179.6 ± 8.8 




4645.79103 


181.6 ± 5.0 


FORS1 


4645.79259 


185.0 ± 3.1 




4692.51274 


28.4 ± 13.4 


TWIN 


4692.52696 


30.7 ± 14.3 




4693.42472 


227.7 ± 11.9 




4693.47199 


137.2 ± 20.6 




4696.49294 


42.0 ± 9.3 






1 nn c a. mo 
l / / .o ± lu.y 




4696.60171 


171.6 ± 9.8 




4979.61251 


90.6 ± 6.8 




4979.65127 


-32.8 ± 21.3 




5035.46989 


207.0 ± 15.3 




5035.49811 


210.0 ± 14.3 




5036.50301 


83.0 ± 15.5 




5037.52184 


-25.0 ± 10.8 




5037.59514 


16.0 ± 15.0 




4758.55029 


206.1 ±22.5 


EFOSC2 


4758.55416 


202.4 ± 25.3 




4758.55803 


212.2 ± 22.5 
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Table B.7. J2256+0656 



mid-HJD 


RV [kms '] 


Instrument 


1 A CA f\f\f\ 

-2450000 






3710.557488 


11.8 ± 10.5 


SDSS 


3710.571464 


1.1 ± 8.6 




3710.588935 


-15.5 ± 9.7 




4694.610760 


78.0 ± 6.0 


TWIN 


4694.676593 


79.3 ± 1.0 




4694.689883 


87.1 ±6.0 




5035.48389 


-140.9 ± 23.8 




5035.50848 


-96.8 ± 15.6 




5036.51782 


95.0 ± 15.0 




5036.61657 


29.9 ± 10.7 




5036.56824 


81.7 ± 11.2 




5037.53430 


-129.6 ± 17.9 




JUj / . J^f / o^f 


nn ft x 1 1 i 

- IZU.o ± 11.1 




5037.55965 


-134.9 ± 11.0 




5037.60833 


-116.5 ± 11.0 




5038.56302 


96.6 ± 10.6 




5038.57506 


115.1 ± 12.7 




5038.62975 


104.1 ± 15.3 




5038.64304 


105.0 ± 11.8 




5038.65573 


104.1 ± 13.1 




5039.64785 


-123.9 ± 12.5 




5048.981437 


16.4 ± 4.2 


GMOS 


5049.002417 


28.5 ± 4.3 




5049.023396 


43.5 ± 4.2 




5049.044376 


51.0 ± 4.3 




5077.83332 


129.0 ± 5.6 




5077.85483 


116.5 ± 5.3 




5068.53076 


-36.5 ± 3.4 


ISIS 


5068.53111 


-30.2 ± 3.5 




5068.53826 


-26.5 ± 2.3 




5068.54815 


-10.3 ± 2.7 




5068.55536 


-12.8 ± 3.2 




5068.56257 


-1.5 ±2.3 




5069.53785 


41.6 ± 3.7 




5069.54507 


48.0 ± 3.8 




5069.55228 


45.8 ± 2.3 




5069.56913 


48.2 ± 3.6 




5069.57635 


26.9 ± 3.8 




5069.58358 


18.2 ±3.8 




5070.62493 


-48.4 ± 4.0 




5070.63214 


-35.8 ±4.0 




5070.63942 


-21.1 ±4.1 




5070.65603 


-20.0 ± 2.6 




5070.66324 


-1.5 ±4.5 




5070.67045 


2.6 ± 4.6 




5071.39970 


8.8 ±4.8 




5071.39997 


21.3 ±4.2 




5071.40718 


33.8 ±3.6 




5071.42392 


50.5 ± 3.8 




5071.43114 


50.8 ±3.9 
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